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Comparing geomagnetic data with data oh tropospheric and stratospheric circulation 
characteristics, we find a statistically highly significant shrinking in areal extent of the strato- 
spheric vortex from the third to the eighth day following a “geomagnetic storm.”. The meridio- 
nality of the 30 640-m contour line at 10 millibars increases markedly from- 5 to 8 days after 
the storm. 

During the contraction of the polar vortex edge, the mean height of the vortex central 
contour decreases only slightly. This indicates that a stratospheric warming event is associated 
with a steepening of the contour gradient rather than a warming over the entire area of the 
stratospheric polar vortex. ' 

The troposphere reacts to these weak, but nevertheless significant, stratospheric warming, 
events by a shrinkage of the area of the 500-miUibar , cold air pool. This shrinkage com- 
mences ,ahout 3 days after the stratospheric warming. 

Our investigation also indicates that the energy input into the stratosphere that is 
received in conjunction with the geomagnetic disturbance has' to come at a propitious time, 

' that is, when the stratospheric-tropospheric circulation system is not already ' undergoing a 
major readjustment because of -ah inherent dynamic instability. It can be shown that the 
observed warming of the stratosphere that follows a geomagnetically disturbed key day can- 
not be explained by simple radiation absorption. 


The complex reaction of. the atmosphere to 
solar geoma^etic activity has become the sub- 
ject of an increasing number of research studies. 
Macdonald, and, Roberts (1960), found that 300- 
millibar troughs that enter or move info the Gulf 
of Alaska ' were amplified several days after 
Earth ‘was bombarded with unusually intense solar 
corpuscular • emission. Macdonald and Roberts 
(1961) and Twitchell (1963) obtained similar 
results of trough intensification at the 500-miUibar 
level.' ' 

Reiter and Macdonald (1973) indicated that 
fluctuations in the area of the tropospheric cold 
pool (r< — 30° C at 500 millibars) and in the 
size of the polar vortex at 10 millibars are coupled 
by a feedback mechanism. They found that sud- 
den warmings in the stratosphere tend to precede 
warmings in the troposphere, and a portion of this 
paper will investigate this stratospheric forcing 


further. Roberts and Olson (1973) indicated that * 
300-milIibar troughs over North America tended 
to intensify with a lag time from a geomagnetic 
event to maximum vorticity development of about 
5 to 7 days. They define a geomagnetic event as 
a daily planetary geomagnetic activity 'index -/Ip, 
greater than or equal to 15 along with an increase 
of this, value over, the previous daily value at 
least as large as the nionthly average value of A^. 

THE GEOMAGNETIC, STRATOSPHERIC, 
AND TROPOSPHERIC DATA AND THEIR 
INTERCOMPARISONS , 

The superposed epoch method was employed 
to investigate a possible relationship between geo- 
magnetic activity and both the wintertime strato- 
spheric polar vortex and the tropospheric cold 
pool. This method compares two sets of data: 
key events are parameterized and selected from 
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one set, and the mean" action or reaction of the 
other set surrounding these key events is noted. 
'In this paper, 29 days surrounding each key event 
are used in each single epoch. These range from 
the 14th day preceding the event to the 14th day 
following it. These dates are noted as D-u, D- 13 , 

. . . , D-i, Do, Di, , Du- The key event occurs 
on Do. 

Specifically, we employed a set of geomagnetic 
activity data to be used in determining the key 
events. We developed two separate sets of data of 
“reacting” events: one dealing with the polar 
troposphere and the other with the polar strato- 
sphere. These three sets of data will be described 
first, and their comparisons and results using the 
superposed epoch method will follow. 

To develop an objective method for determin- 
ing a sudden increase in geomagnetic activity, we 
used the daily planetary geomagnetic activity 
index Ap, as published by the National Geophysi- 
cal and Solar-Terrestrial Data Center. This is a 
daily global index of geomagnetic activity and is 
generally considered to be linear to its severity.' 
Key dates of this activity, called “geomagnetic key 
dates,” were selected according to two criteria: 
The daily Ap value must be greater than or equal 
to 15, and the increase from the previous daily 
value must be at least as large as the monthly 
average value of Ap. These are the same two 
criteria used in the paper by Roberts and Olson 
(1973). The key dates cover 17 yr from 1953 
through 1969 and therefore are available for all 
winters for which we have tropospheric and strat- 
ospheric data available. 

Our set of data for the stratosphere parame- 
terizes the size and convolution of the polar vortex 
at 10 millibars. It is identical to that used in the 
previous study by Reiter and Macdonald (1973). 
The 30 640-m contour at this pressure level gen- 
erally lies near the edge of the polar vortex during 
the months from November through March. The 
latitude value of this contour at 30° longitude 
intervals is noted for each day, giving 12 such 
values. The mean of these latitudes gives a rough 
idea of the daily areal extent, although not of the 
intensity, of the vortex. The standard deviation of 
these values gives an indication of the convolution 
or ellipticity of the vortex. For each day in the 
12 cold seasons (November through Match) 


1957-58 through 1968-69, we obtained a mean 
latitude value as well as a standard deviation 
value for' this contour line. 

The tropospheric data deal with the daily size 
of the 500-millibar cold pool. Generally, the 
—30° C isotherm lies near the polar front at this 
level, and the area enclosed by this isotherm 
should give an indication of 'the areal extent of 
the cold pool. We planimetered the area enclosed 
by, this isotherm from maps published by the U.S. 
National Weather Service for each day in 10 cold 
seasons, 1953-54 through 1962-63. Values for 
two of the seasons, 1961-62 and' 1962-63, were 
taken from operational charts while the others 
were taken from the Daily Series Synoptic Weath- 
er Maps published by the U.S. National Weather 
Service. Portions of this area that occasionally 
broke away from the main cold pool were disre- 
garded unless they “rejoined” the pool at a later 
time. This data set consists of the daily area of the 
500-millibar cold pool in arbitrary units., 

Comparisons of Geomagnetic Data With 
Stratospheric and Tropospheric Data 

First let us compare the geomagnetic key dates 
with the mean latitude and standard deviation of 
the polar vortex, our stratospheric data. Ninety- 
eight key dates were selected from nine cold sea- 
sons, 1960-61 through 1968-69. The mean 
values of these two sets of stratospheric data for 
the 98 epochs surrounding the key events are 
shown in figure 1. Note the significant increase 
in mean latitude of the 30 640-m contour, indi- 
cating a shrinkage of the polar vortex, from the 
third to the eighth day following the geomagnetic 
event. The Wilcoxon Rank Sum Test shows that 
the Di through Du mean latitudes are statistically 
separate from the D-u through D-i means at the 
99-percent significance level. Most perplexing is 
the slight increase in mean latitude along with a 
corresponding sharp increase in standard deviation 
preceding the key date. To investigate this situa- 
tion, we reduced our key dates to only those 
which were preceded by at least nine nonkey 
dates. This eliminates the “preevent” compound- 
ing effects of sequences of key events. Forty key 
dates met this new criterion, and the mean values 
of the mean latitude and standard deviation for 
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Figure 1. — Superposed epoch averages of the daily 
mean latitude "h (top) and the daily standard devia- 
tion a (bottom) of the 30 640-m contour line at 10 
millibars surrounding key geomagnetic dates. Data 
averaged were taken from 98 cases in 9 cold seasons 
(November through March) for the years 1960-61 
through 1968-69. 

these epochs are show in figure 2. It was noticed, 
however, that a sudden breakup of the polar 
vortex circulation occurred during two of these 
epochs: the mean latitude of the 30 640-m con- 
tour fluctuated by as much as 20° in one day in 
these two cases. The mean latitudes of these two 
individual epochs are shown in figure 3. After 
eliminating these sequences, we are left with the 
mean values of 38 epochs, which are shown in 
figure 4. Note the rapid increase in mean latitude 
from D 3 through D?. Also, the standard deviation 
of the vortex jumps most markedly from D 5 
through Dg. These figures indicate that a 4- to 
5-day shrinkage of the polar vortex follows a 
key geomagnetic date by about 3 days, with a 
slight increase in the ellipticity of, or meridional 
transport by, the polar vortex later in the period 
of the shrinkage. 

Returning to the 98 original epochs and taking 
them individually, we tried to determine the statis- 
tical significance of the Dr through Dn mean lati-^ 
tudes compared with some prekey event values. 



Figure 2 . — Key geomagnetic dates that were preceded 
by at least nine nonkey geomagnetic dates (40 cases) 
in 9 cold seasons (November through March) for 
the years 1960-61 through 1968-69. 

Specifically, we used the D -10 through D-i mean 
latitudes for the preevent data, giving a total of 
15 values to be compared for each epoch. A 
simple rank sum test was used to compare these 
two sets of data and to determine the statistical 
significance of their separation. In 52 of the 98 
epochs, the mean latitude of the Dr through Dn 
data is greater than the preevent values at the 95 
percent significance level. In other words, in more 
than half of the key epochs, this Dr through Du 
increase in mean latitude following the key event 
is significant. 

Three seasons with stratospheric and geomag- 
netic data (1957-58 through 1959-60) remain, 
and we used these data to determine whether the 
same, trend will develop from new independent 
data. Thirty-one key geomagnetic dates were 
chosen from this sample, and the results of the 
superposed epoch method of mean latitude and 
standard deviation determination are shown in 
figure 5. Again we selected only those key dates 
that were preceded by at least 9 nonkey dates, of 
which there were 14, and the results of the super- 
posed epochs for these events are shown in figure 
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Figure 3. — The daily mean latitude values of the 
30 640-m contour at 10 millibars surrounding the key 
geomagnetic dates of January 30, 1963 (solid line), 
and February 10, 1973 (dashed line). 

6. Note a similar trend toward an increase in 
mean latitude following the geomagnetic event 
(in this case from 6 to 8 days following the key 
date). The large increase in standard deviation 
preceding the key date is due mostly to a single 
event, while the increase preceding Ds is more 
general. 

We also tried to determine a mean 500-millibar 
cold pool response surrounding similar geomag- 
netic events. Because the tropospheric data and 
the stratospheric data cover different seasons, the 
key dates are not exactly the same; however, the 
criteria used in selecting them remain identical. 
The 10 cold seasons that were used ran from 
1953-54 through 1962-63, and 113 days were 
selected as key geomagnetic dates from this 



Figure 4. — The daily mean latitude values at 10 millibars 
(38 cases). 



Figure 5. — Superposed epoch averages of the daily mean 
latitude 4> and the daily standard deviation <r of the 
30 640-m contour line at 10 millibars surrounding key 
geomagnetic dates for the 1957-58 through 1959-60 
cold seasons (31 epochs). 
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Figure 6. — Key geomagnetic dates that are preceded by 
at least nine nonkey dates for the 1957-58 through 
1959-60 cold seasons (14 epochs). 


period. The mean values of the area within the 
—30° C isotherm surrounding the key dates are 
shown in figure 7(a). No statistically significant 
variation can be determined from these data. 

I . ' • 

Selecting only those key dates that were preceded 
by at least nine nonkey dates, we noted the mean 
area variations that are given in figure 7(6). 
Again, no significant variation is apparent. 

Sector Events 

Occasionally, and often at the time of a geo- 
magnetic event, the orientation of the interplan- 
etary magnetic field switches. Wilcox et al. (1973) 
observed a vorticity minimum in the troposphere 
and lower stratosphere north of 20° N latitude 
about 1 day following the passage of a sector 
boundary. No overlap of our tropospheric and 
sector data was available, but we wanted to deter- 
mine whether such a switch had an effect on the 
stratospheric polar vortex at 10 millibars. Forty- 
two dates of this switch, whether from positive to 
negative or. vice versa, were selected from the cold 
seasons 1963-64 through 1968-69. These were 
called sector key dates, and the superposed epoch 
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Figure 7. (a). Superposed epoch averages of the daily 
area (in arbitrary units) of the cold poo! (T^— 30° C) 
at 500 millibars surrounding' key geomagnetic dates. 
Such dates (113 in all) were selected from November 
through March in the seasons 1953-54 through 1962- 
63. (fc). Superposed epoch averages of the daily area 
(in arbitrary units) of the cold pool (T^~30° C) at 
500 millibars surrounding key geomagnetic dates. Key 
dates include only those preceded by at least nine non- 
key dates (45 cases) and were selected from Novem- 
ber through March in the seasons 1953-54 through 
1962-63. 

method was used to determine a mean strato- 
spheric reaction surrounding these dates. The 
mean of the 30 640-m contour mean latitude and 
the mean of its standard deviation surrounding 
these key events are shown in figure 8 . 

Note the slight decrease in mean latitude 
(expansion of the polar vortex) following the key 
date, with relatively lower values from D 3 through 
Z>T. Using a simple rank sum , test, we compared 
the values for these 5 days with those of the D -10 
through D-i segment separately for each of the 42 
sequences. In 14 of the cases, the D 3 through Dr 
sample was lower than the prekey date sample at 
the 95-percent significance .level. In 16 of the 
cases, however, this Ds through Dr sample was 
actually greater than the prekey date sample above 
the 95 percent significance level. Thus we could 
establish no statistically significant trend. 
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Figure 8.^ — The superposed epoch averages of the daily 
mean latitude d’ and the daily standard deviation a of 
the 30 640-m contour at 10 millibars surrounding 
sector key dates. Forty-two cases were included from 
November through March in the seasons 1963-64 
through 1968-69. 

Tropospheric Response to the Stratosphere 

We have shown that there appears to be a 
stratospheric reaction to geomagnetic activity, but 
there appears to be no similar significant response 
in the troposphere. Reiter and Macdonald (1973) 
indicated that the troposphere reacts to sudden, 
strong warmings in the stratosphere and that these 
tropospheric warmings tend to occur about 2 days 
later. (See figure 9.) We wanted to include the 
effects of weaker and less sudden warmings in 
the stratosphere in this study, however. Using our 
stratospheric data for the six seasons in which it 
overlapped the tropospheric data (1957-58 
through 1962—63), we took every possible 9-day 
sequence in each season and separated it into 
three 3-day sequences. Key stratospheric warm- 
ing events were determined in the following man- 
ner: the 30 640-m contour mean latitude in the 
second 3-day sequence must be greater than the 
mean of the first 3-day sequence by 2° of latitude 
or more, and similarly the mean of the third 3-day 
sequence must also be greater than the second by 
2° or more. Key dates were arbitrarily called the 
fifth day (the middle day) of the 9-day sequence, 
and 52 such sequences in the six seasons met both 
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Figure 9. — Superposed epoch averages of four cases of 
stratospheric vortex breakdown measured by an in- 
crease in the mean latitude of the 30 640-m contour 
at 10 millibars (/op) and the mean area (arbitrary 
■ units) of the cold air (7^-30° C) at 500 millibars 
(i>o//om).- (From Reiter and Macdonald, 1973.) 

criteria. Using the superposed epoch method, we 
determined the mean response of the tropospheric 
cold pool area surrounding these key dates. The 
mean values of the polar vortex mean latitude 
(the controlled event) and the 500-millibar cold 
pool area are given in figure 10. Note the shrink- 
age of the cold pool following the stratospheric 
warming, with the most significant shrinkage 
beginning about 3 days after the stratospheric 
warming. To test the statistical significance of this 
decrease in area, we again used a simple rank sum 
test separately for each of the 52 sequences. We 
compared the area values of the D-s through Z)-i 
sequence with those of the Dg through D 12 se- 
quence. In 32 of the 52 epochs, the latter sample 
was statistically less than the former sample at 
the 95-percent significance level or better. In 40 
of the cases, the numerical mean of the Dg through 
Di 2 sequence was less than the mean of the ear- 
lier sequence. This confirms a forcing upon the 
tropospheric cold pool size by stratospheric warm- 
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Figure 10. — Superposed epoch averages of the 30 640-m 
contour mean latitude ^ at 10 millibars surrounding 
an increase in mean latitude of 4° or more in 9 days 
(lap), and supeiposed epoch averages of the area of 
the cold air (T<— 30° C) in arbitrary units surround- 
ing such events (bottom). 


ings that are weaker than those discussed by 
Reiter and Macdonald (1973). ^ 

We speculate that the reason that no tropo- 
spheric response to geomagnetic activity could be 
shown directly is that the intermediary action of 
the stratosphere tends to mask this effect over the 
time scales considered here. This would cause the 
tropospheric reaction to be spread over a greater 
length of time with respect to the key geomagnetic 
date; therefore, it would be more dfficult to detect 
in a statistical sense. 

The results presented in this section indicate 
that the stratosphere responds more significantly 
to geomagnetic activity than does the troposphere, 
and that the resulting stratospheric warming is in 
turn forced upon the troposphere. This forcing has 
been the subject of. several earlier papers (Austin 


and Krawitz, 1956; Reiter and Macdonald, 1973; 
Teweles, 1958). ‘ ^ ^ . 

POSSIBLE MECHANISMS 
Polar Vortex Center 

Before determining the mechanism that brings 
about the shrinkage of the polar vortex discussed 
in the preceding section, it is important to examine 
the fluctuations of the vortex center surrounding 
such warming events. If the center contour at 10 
millibars shows a marked increase at the time that 
the edge of the vortex shrinks, a mechanism of 
large-scale subsidence would suggest itself. A 
schematic indication of a typical event of this type, 
if it exists, is shown in figure 11. On the other 
hand, if the center contour remained essentially 
at the same value or became numerically less 
during shrinkage, a steepening of the contour gra- 
dient near the edge of the vortex would be asso- 
ciated with a contraction of the vortex edge. Some 
mechanism such as mass importation or warming 
only along a rather narrow belt would be indi- 
cated. Figure 12 shows a schematic interpretation 
of an event of this type. 

We examined the fluctuations in central con- 
tour value during a 29-day epoch surrounding a 



Figure 11. — Meridional cross section of the 10-millibar 
surface surrounding an increase in mean latitude 
(shrinkage of the polar vortex) of the 30 640-m con- 
tour, if it is associated with large-scale warming or 
subsidence. The solid line represents the 10-miIIibar 
heights preceding the shrinkage, and the dashed line 
represents height values following the shrinkage. 
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Figure 12,— Meridional cross section of’the ip-millibar 
surface surrounding an increase iri mean latitude 
(shrinkage of the polar vortex) of the 30 640-m con- 
tour, if it is associated with a steepening of the contour 
gradient along the vortex edge. The solid line repre- 
sents the 10-millibar heights preceding the shrinkage, 
and the dashed line represents height values following 
the . shrinkage, ■ 

contraction, of the vortex edge. As before, we used 
the criterion in which the mean latitude of the 
30 640-m contour at 10 millibars increased by 4° 
or moire in 9 days using the method with the 3- 
day means described earlier. The superposed 
epoch method was employed with the key date 
chosen again to be the middle day of such 9-day 
sequences. In the 12 seasons for which we have 
10-miIlibar data, seventy-six 9-day sequences met 
the criterion. The means of the 30 640-m mean 
latitude values for these events are shown in 
figure 13. The means of the central contour value 
at 10 millibars during these epochs are also shown 
in figure 13. Note that no increase in height of 
this pressure surface, is even remotely suggested; 
in fact, a mean decrease of about 20 m is. implied. 
On the basis' of these results, we can rule out any 
mechanism that' promotes large-scale subsidence 
as being responsible for a shrinkage of the polar 
vortex. We are forced to rely on a mechanism 
that causes a steepening of the contour gradient 
(on a constant pressure surface) near the edge of 
the polar vortex to bring about the observed 
contraction. • , ■ , - 

• One .possibility of warming the polar vortex 
edge at 10 millibars would be through collisional 




Figure ! 3. — Superposed epoch averages of the 30 640-m 
contour mean latitude at 10 millibars surrounding 
an increase in mean latitude of 4° or more in 9 days 
{top) and superposed. epoch averages of the value (in 
meters) of the polar vortex central contour at 10 
millibars {bottom). 

excitation and ionization of the atmospheric mole- 
cules during the geomagnetic storm; i.e., throng 
direct, absorption of energy. Certainly the fact 
that auroras occur along a latitude belt near the 
polar vortex edge gives impetus -to an investiga- 
tion of this possibility. We' will present some calcu- 
lations showing that this mechanism cannot supply 
the required energy to bring, about the observed 
contraction. ' • , , 

According to Matsushita and Campbell (1967), 
we can assume that the auroral absorption takes 
place primarily in a latitude band 10° wide, 
averaging, 5000 km in length in both hemispheres. 
The rate of dissipation resulting • from' auroral 
processes during a magnetic storm is about 10^’ 
to 10'®- erg • s"'. The area of one of these, bands 
is about 5.6 X 10'® cm^,and we will assume that 
10'® erg' • s~' are absorbed over one of these 
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bands during a magnetic storm. A cursory exam- 
ination of the contour gradient at 10 millibars 
near the polar vortex edge in midwinter yields a 
mean contour gradient of about — 80 m per de- 
gree latitude, shown schematically in figure 14. 
If we assume uniform -heating of a 10° latitude 
band (from 50°. to '60° N) only, a 4° increase in 
mean latitude of the 30 640-m contour line would 
require a uhiform 320-m increase in height of 
the 10-millibar surface over this latitude band. If 
this increase is due. totally to heating in the 30- 
to 10-millibar layer, the calculations shown in 
appendix A indicate a required mean warming of 
about 10° C in this layer. Also in appendix A, 
calculations of the energy required to carry on 
this heating compared with the energy available 
from a long (10^ s) geomagnetic event show that 
simple absorption and redistribution of the auroral 
energy could not possibly account for the noted 
heating. 

DISCUSSION 

■ It is apparent that simple abso^tion of the 
radiative energy associated with a geomagnetic 
storm cannot account for the observed warming 
af 10 millibars following such an event. Some 
mechanism involving the dynamics and transport 
processes along the vortex edge should.be investi- 



Figure 14. — A schematic diagram of 10-millibar sur- 
faces with latitudinal gradient of — 80 m per degree 
latitude. ' ' ' . 


gated. In particular, adiabatic sinking motion and 
eddy transport processes in the area might account 
for the observed warming. Calculation of the 
adiabatic subsidence in the 30- to 10-millibar 
layer required to produce such a- warming are 
shown in appendix B. The result (0.14 cm ^ s~^) 
is within the realm of variability in vertical 
motion at 50 millibars reported by Mahlman 
(1966). He indicates that mean vertical motion 
during a. “stratospheric warming” changed from 
—0.06 cm • s"^. preceding the period to —0.14 
cm • s“' during it. The increase in standard devia- 
tion of the 30 64Q-m contour at 10 millibars (see 
fig. 4) indicates that the. effect of eddy transport 
processes is increasing after a geomagnetic key 
date, and this too may- account for some of the 

observed warming. 

• r ‘ 
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APPENDIX A— CALCULATIONS OF 
ENERGY REQUIRED FOR 
STRATOSPHERIC WARMING VERSUS 
AURORAL ENERGY 

(1) Assume a mean temperature of 218 K 
(—55° C) in the 30- to 10-millibar layer. 

(2) Given the formula from the Smithsonian 
tables: 

= 67.442(273.16 + log^ 

where 

A$ = thickness of the layer, geopotential meters 
Tgpm) ,1 . 

t'mv = mean adjusted virtual temperature of the 
layer, °C . 

Pi = pressure at the base of the layer 
Pt = pressure at the top of the layer. 

(3) Using this formula with the values given 

in(l), . . . 

A$ = 7020 gpm 

(4) If we increase the thickness of this, layer 
by 320 gpm and reapply the equation in'(2), 

fmv = -45° C • ■ ■ • 
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(5) Therefore, corresponding to an increase of 
320 gpm in the 30- to 10-millibar layer, the mean 
virtual temperature must increase by 10° C. 

(6) From the text, we had assumed that the 
'area of the latitude band in which auroral energy 
is absorbed is 5.6 X 10^® cm^. 

(7) The mass of air in the 30^ to 10-millibar 
layer over this band is 

(20 g . cm-^)(5.6 X 10“ cin^) = 1.1 X 10“ g 

(8) The specific heat of air Cp is given as 

10 erg • gr^ • K'^ 

(9) The energy required to bring about this 
observed warming is equal to the total mass to be 
heated multiplied by the specific heat of the mass 
multiplied by the change in temperature required, 
from (7), (8), and (5): 

Energy required = 

(1.1 X 10“ g)(10® erg vg"- ■ K-^(10 K) 

= 1.1 X 10“ erg 

(10) From Matsushita and Campbell (1967), 
assume that the energy of an auroral absorption is 
10“ erg • s~^. 

(11) Assume that this strong absorption lasts 

3 hr or 10* sec. 

(12) Then the total energy involved in the 
aurora is 

( 10“ erg . s-*)(10* s) = 10“ erg 

(13) Comparing the results from (9) and (12), 
note that the energy involved in an aurora is 
much less than is required to produce the noted 
heating. 

APPENDIX B— CALCULATIONS OF 
SUBSIDENCE REQUIRED FOR 
STRATOSPHERIC WARMING 

Assume a 4° increase in mean latitude of the 
,30 640-m contour at 10 millibars and assume that 
this is brought about by the 10 K warming in the 
30- to 10-millibar layer noted in appendix A. 

Differentiating Poisson’s equation and holding 
de = 0 where P = 20 millibars and T = 223 K, 
let dr = + 10 K: 

de = dr(i^^ - xr (1000)'" dP 

dP = 3.1 millibars 


Using the hydrostatic approximatioii, this cor- 
responds to a change of about 1070 gpm. 

Therefore a parcel of air that sinks adiabatically 
from the 20-millibar level at T = 223 K and 
warms 10 K must experience a change in geopo- 
tential of ~1070 gpm. 

If this change in geopotential is experienced 
over a period of 9 days (7.78 X 10® s), then the 
mean vertical motion that accounts for this warm- 
ing is about —0.14 cm • s"'. 
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DISCUSSION 

SHAPIRO: Could you define a little more precisely 
the nature of your, magnetic key day selection? 

MACDONALD: We used a planetary Ap index to 
determine these key dates. It had to be at least 15, and 
the increase over the previous day had to be greater 
than or equal to, the mean monthly Ap value. 

SHAPIRO: That is similar to what Roberts has done. 
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MACDONALD: That’s exactly the same criterion he 
used, yes. . _ . - . 

AJD^N: Have you made ahy analysis on whether the 
polar vortex ever breaks up in association with geomag- 
netic activity? 

MACDONALD: Yes; in fact, it did break up. A 
breakup occurred near a key date twice, I believe. We 
excluded such data from these charts to avoid the mask- 
ing of any other values that we observed from, say, the 


other 38 key dates; but we only had 12 yr of these data 
and. we. could detect no real correlation, with, for exam- 
ple, a massive breakup of the polar vortex following 
that key date. 

QUESTION: What time of the year did the breakup 
occur? 

MACDONALD: There were two breakups that oc- 
curred near a key date, and they were both in January. 
Our data run from November through March. 



